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1. Introduction 

The concept of pointer-based simultaneous measurements of conjugate observables is an 
indirect measurement model, which allows to dynamically describe the properties of a 
simultaneous quantum mechanical measurement process. Additional to the system to be 
measured (hereafter just called system ), the model introduces two additional systems 
called pointers, which are coupled to the system and act as commuting meters from 
which the initial system observables can be read out after a certain interaction time. 
In this sense, the pointers represent the measurement devices used to simultaneously 
determine the system observables. Pointer-based simultaneous measurements date back 
to Arthurs and Kelly [1] and are based on von Neumann’s idea of indirect observation 
[ 2 ]- 

In principle, any pair of conjugate observables like position and momentum 
or quadratures of the electromagnetic field [3], whose commutator is well-defined 
and proportional to the identity operator, can straightforwardly be measured within 
the scope of pointer-based simultaneous measurements. We limit ourselves to the 
measurement of position and momentum in the following. An open pointer-based 
simultaneous measurement [4] also takes environmental effects into consideration by 
utilizing an environmental heat bath in the sense of the Caldeira-Leggett model 
[5, 6, 7, 8], which leads to a quantum Brownian motion [9, 10, 11, 12] of the 
system and the pointers, whereas a closed pointer-based simultaneous measurement 
[1, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22] does not involve any environmental effects. A 
schematic open pointer-based simultaneous measurement procedure is shown in Fig. 1. 

In this contribution we calculate the information entropy of an open pointer-based 
simultaneous measurement and discuss its properties as a measurement uncertainty. 
In particular, we make use of recent results [22, 4], which we extend and generalize. 
In Sec. 2, we present the formal dynamics of open pointer-based simultaneous 
measurements and then use these results to discuss the entropic uncertainty in Sec. 3. 
In the end, we arrive at a generic lower bound of this entropic uncertainty. Note that 
we solely use rescaled dimensionless variables [4] so that h — 1. 
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infer initial system observables 



Figure 1 . Principles of an open pointer-based simultaneous measurement of two 
conjugate observables, e. g., the simultaneous measurement of position and momentum. 
The measurement apparatus consists of two quantum mechanical systems, called 
pointers, which are bilinearly coupled to the quantum mechanical system to be 
measured. Additionally, an environmental heat bath in the sense of the Caldeira- 
Leggett model can disturb both the system and the pointers. After the interaction 
process, one observable of each pointer is directly measured (e. g., the position of each 
pointer) while the system itself is not subject to any direct measurement. However, 
from these measurement results, information about the initial system observables can 
then be inferred. In other words, the final pointer observables act as commuting 
meters from which the initial non-commuting system observables can be simultaneously 
read out. The price to be paid for this simultaneity comes in form of fundamental 
noise terms, which affect the inferred values. The corresponding uncertainties can be 
described by information entropies, which are bound from below. 


2. Open pointer-based simultaneous measurements 

As indicated in the introduction, our model of open pointer-based simultaneous 
measurements consists of a system particle to be measured with mass Mg, position 
observable Xg and momentum observable Pg, which is coupled bilinearly to two pointer 
particles with masses Mi and M 2 , position observables X\ and X 2 , and momentum 
observables and P 2 , respectively. Both the system and the pointers are bilinearly 
coupled to an environmental heat bath, which consists of a collection of N harmonic 
oscillators with masses position observables < 71 ,..., <?jv, an d momentum 

observables ki,... ,kN- In this section, we first present the general Hamiltonian for this 
model and then briefly discuss the resulting dynamics. 

2.1. Hamiltonian 

The general Hamiltonian for our model reads 


= Pf r ee + H int (t) + Pbath(^) 


( 1 ) 


and therefore consists of three parts. First, the free evolution Hamiltonian 
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which simply describes the dynamics of the undisturbed system and pointers. Second, 
the interaction Hamiltonian 


H int (t) = C s (t)Xi + t)Xl + C 2 {t)Xl + (X s , Ps)C(f)(X 1? x 2 , A, A) T , (3) 

which describes possible quadratic potentials with the coupling strengths Cs(t), Ci(t), 
and C 2 (t), respectively, as well as bilinear interactions between the system observables 
and the pointer observables via the 2x4 coupling matrix C (t). These interactions are 
necessary for an information transfer between system and pointers and are therefore a 
prerequisite of pointer-based simultaneous measurements. The existence of quadratic 
potentials is on the other hand not essential, but may be reasonable from a physical 
point of view when regarding confined particles. One possible interaction Hamiltonian 
would be the interaction Hamiltonian of the classic Arthurs and Kelly model [1], which 
can be written as H int = k(X s A + AA) with an arbitrary coupling strength n ^ 0. 

Lastly, Eq. (1) contains the bath Hamiltonian [5, 6, 7, 8] 

Abath(f) = ^k T m -1 k + ^q T cq + q T g(t)(X s , X u X 2 ) T , (4) 

which describes the independent dynamics of the bath particles with the Nx N diagonal 
mass matrix m containing m 1; ..., m,v, and the N x N symmetric and positive definite 
bath-internal coupling matrix c; as well as the coupling of system and pointer positions 
to the bath positions with the coupling strength 


s{t) = g(t) g, (5) 

which consists of the time-dependent scalar g(t) and the time-independent N x 3 
matrix g. To simplify our notation, we make use of the vectorial bath positions 
q T = (g lf ..., q N ) and the vectorial bath momenta k 7 = (A, ■ ■ ■, k N ). 

In particular, since all of the following calculations only rely on the bilinear structure 
of the Hamiltonian, Eq. (1), we do not further specify the coupling strengths and 
therefore consider quite general measurement configurations. Note that the time- 
dependencies of the coupling strengths in Eqs. (3) and (4) allow us to design specific 
coupling pulses for system-pointer interactions [20] or switch-on functions for the bath 
[10, 23], 

2.2. Dynamics 

To determine the complete system and pointer dynamics, it is necessary to solve the 
coupled Heisenberg equations of motion 
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( 6 ) 


which take on the form of a Volterra integro-differential equation [24] after explicitly 
solving the dynamics of the bath particles [25, 10, 11]. The general solution of Eq. (6) 
can formally be expressed by means of a resolvent [26, 27] and leads to system and 
pointer positions and momenta, which are linearly propagated from their initial values 
(i.e., the homogeneous solution) and are affected by an additive noise term (i. e., the 
inhomogeneous solution). Explicit analytical solutions, which are naturally of the 
same structure as the formal solution, can only be found for specific choices of the 
Hamiltonian, Eq. (1). Otherwise, numerical approaches [28, 29] might be necessary. 
Interestingly, since we strive after a general discussion of the dynamics, the ensured 
existence of a solution and its formal structure is sufficient for all further considerations. 
We assume in this context that possible unphysical artifacts of the modeling [30, 10], 
e. g., renormalization problems, have been treated adequately [4], Without loss of 
generality, we choose t — 0 as the initial time. 

For the description of a pointer-based simultaneous measurement, we do in fact not 
need knowledge about the complete system and pointer dynamics. As outlined in Fig. 1, 
the measurement process provides that the two pointers are being measured after an 
interaction time t in such a way, that we either read out the position or the momentum 
of each pointer. Consequently, we have knowledge about either X\ (t) or A (t) of the 
first pointer, and, additionally, about either X 2 (t) or P 2 (t) of the second pointer, which 
is a total of four different possible measurement combinations. To summarize these 
measurement combinations, we define the measurement vector w (t), which consists of 
the two pointer observables chosen to be read out, e.g., w (t) = (Xi(t), X 2 (t)) T when 
measuring both pointer positions. As also outlined in Fig. 1, the information gained from 
measuring the observables in the measurement vector allows us to infer the initial system 
observables A's(O) and A(0). For this reason, the connection between the measured 
observables after a certain interaction time t and the initial system observables builds 
the framework for a description of pointer-based simultaneous measurements. 

With this purpose in mind, we can extract 




from any (formal) solution of the Heisenberg equations of motion, Eq. (6). Here we have 
introduced the so-called generalized inferred observables 



( 8 ) 
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which are given by a rescaled measurement vector w (t) and can on the other hand be 
understood as the effectively measured observables from which the system observables 
can be directly read out. Since the two measured pointer observables in w (t) initially 
commute and are subject to a unitary evolution, one has 


[m,nt)\ = o, (9) 

which means that the inferred observables can be determined simultaneously. 

Furthermore, Eqs. (7) and (8) contain the coefficient matrices A (t), B(t), and 
Aft, s ) with 0 < s < t, for which we do not need to give an explicit expression. 
In general, they can be straightforwardly calculated from the resolvent of the formal 
solution of the complete system and pointer dynamics. An exemplary calculation can 
be found in Ref. [4]. We also make use of the initial value vector 


j = (x 1 (o),x 2 (o),A(o),A(o)) t , 


and the stochastic force [30] 


( 10 ) 


£(t) 
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with the symmetric bath frequency matrix [11] 


( 11 ) 


U) 



cm 


1 
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( 12 ) 


In particular, the stochastic force results from the homogeneous bath dynamics and 
describes the noisy influence of the bath on the measurement results, Eq. (4). The 
symbol * in Eq. (7) represents the integral 


X 

{f* 9 )(x) = Jdyf(x,y)g(y) (13) 

0 

for two arbitrary functions f{x,y ) and g{x). I 11 case of f(x,y ) = f(x — y), Eq. (13) is 
called a Laplace convolution. 

The central aspect of the dynamics contained in Eq. (7) can be understood when 
considering the expectation value 




( 14 ) 
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of Eq. (7), where the brackets refer to the mean value with respect to the initial state. 
Most importantly, Eq. (14) shows that the knowledge about the first moments of the 
inferred observables, Eq. (8), which can be determined by measuring the two pointer 
observables contained in the measurement vector w (t), allows us to infer the first 
moments of the initial system observables. The second and third term on the right- 
hand side of Eq. (7) simply add the shift 

s(f) = B(t)(J> +(A* <£»(*) (15) 

to this relation. Assuming a suitably separable initial state, this shift is determined by 
the initial state of the pointers and the bath, but independent of the state of the system. 
In this sense, our model is only of any practical purpose if the initial pointer and bath 
states are sufficiently well-known to determine s(f). Generally, the initial pointer states 
are at our disposal and can therefore be chosen in such a way that the Erst term on 
the right-hand side of Eq. (15) vanishes. It is furthermore reasonable to suppose that 
an environmental heat bath in thermal equilibrium leads to a vanishing expectation 
value of the stochastic force £(£), Eq. (11). Consequently, presuming a vanishing shift 
s(f) seems appropriate for a typical measurement configuration. We will confirm this 
presumption further below for a specifically chosen initial state. Note that the linearity 
of Eqs. (7) and (8), which is essential for the inference process, Eq. (14), is a direct 
result of the bilinear structure of the Hamiltonian, Eq. (1). 

We assume here the non-pathological case that the interaction Hamiltonian H int (t), 
Eq. (3), and the measurement vector w(f) are chosen in such a way that inferring system 
observables from the pointers is possible in the Erst place, which is defined by the 
existence of the coefficient matrix A (t), Eq. (8). In other words, we require a sufficient 
information transfer from the system to the pointers. For example, for the classic 
Arthurs and Kelly model [1] with the interaction Hamiltonian mentioned in Sec. 2.1 
and no environmental heat bath, measuring both pointer positions leads to an existing 
matrix A(f) for t > 0, but measuring both pointer momenta does not [20]. 

Seeing now the role played by the inferred observables, we can quantify the 
uncertainty of a simultaneous pointer-based measurement with the help of the so-called 
noise operators [31] 


( N x (t)\ (X(t) - A' s (0)\ 

yvp(i)) \v(t) - P s (o)) 


B(f)j +(A *«)((), 


(16) 


which are defined as the difference between the inferred observables X{t) and V(t), 
Eq. (8), and the respective initial system observables X s (0) and Pg(0). The second 
equal sign in Eq. (16) directly follows from Eq. (7). The expectation values 



= S W 


N x (t) 
Nvit ) 


(17) 
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of the noise operators correspond to the shift s(t), Eq. (15). As we have mentioned 
above, this shift can typically be presumed to vanish. Particularly interesting for our 
considerations is the symmetrized covariance matrix 


( (iVjM) - (N x (t)) 
\{NMt)) - (Mt)) (M*)) 


(NMt)) - (#*(<)) (jVp(t))\ 

(jv?w) - ( mt )) 2 ) 


of the noise operators, where 


d 2 x {t) 8 xv {t)\ 

8xv(t) J 

(18) 


N xv (t) = ^ (N x (t)N v (t) + N v (t)N x (t)) (19) 

stands for the symmetrized noise operator. To simplify our notation, we have also 
introduced the noise terms 5 x (t) and 8f,(t), which represent the diagonal elements of the 
covariance matrix, and the correlation term 8 x -p(t), which represents the off-diagonal 
elements. Since a precise knowledge of s (t), Eqs. (15) and (17), is necessary for the 
success of the inference process, Eq. (14), the associated variances given by the noise 
terms 8 x (t) and 8%,(t) can be considered as an intuitive measure for the respective 
measurement uncertainty. Indeed, in the scope of variance-based uncertainty relations, 
the noise terms play an important role; see, e.g., Ref. [32] and references therein. Note 
that the noise terms are also referred to as “errors of retrodiction” [16, 17, 18]. We will 
revisit them further below, where we will see that they arise naturally in our entropic 
description of the measurement uncertainty. 

The formal dynamics of the inferred observables X(t) and V(t), Eq. (7), build 
the framework for all of the following considerations. Most importantly, the intimate 
connection between X (t) and V(t), which result from the measured pointer observables 
in w(i), Eq. (8), and the initial system observables Ag(0) and -Pg(O) can clearly be 
seen from Eq. (7). It is this connection which lies at the heart of the pointer-based 
measurement scheme: information about the non-commuting system observables can 
be gathered from knowledge about the commuting inferred observables. At this point it 
seems natural to ask: What is the uncertainty of such an indirect measurement process? 
In the next section, we answer this question with the help of information entropy. 


3. Entropy 

The uncertainty principle [33, 34, 19] is of central importance for quantum mechanical 
measurements. It manifests itself in the form of uncertainty relations [35], usually 
written in terms of variances [36, 37] or information entropies [38, 39]. The concept of 
information entropies goes back to Ref. [40], whereas its general usage in the scope of 
uncertainty relations has been pioneered by Ref. [41, 42, 43]. In the context of closed 
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pointer-based simultaneous measurements [1, 14], the uncertainty principle leads to well- 
known variance-based uncertainty relations [16, 17, 18], which have extensively been 
discussed, e. g., in the scope of phase-space measurements [13] or energy and timing 
considerations [20, 21]. A respective information entropic form has been derived in 
Ref. [15] and improved in Ref. [22], Variances can also be used to describe an uncertainty 
relation for open pointer-based simultaneous measurements [4] . However, in comparison 
with information entropies, variances suffer from two major drawbacks [44, 15, 45]: 
First, they can become divergent for specific probability distributions and second, they 
may not reflect what one would intuitively consider as the “width” of a probability 
distribution. On the other hand, a disadvantage of information entropies is that they 
are usually much more difficult to calculate than variances. This, however, is only a 
technical limitation. Therefore, we concentrate on information entropic uncertainties in 
the following. 

In this section, we first introduce the so-called collective entropy as a total 
measure of uncertainty in the context of open pointer-based simultaneous measurements. 
Choosing a separable initial state then allows us to calculate the marginal probability 
distributions this collective entropy is based on. As a result, we can discuss the structural 
properties of the collective entropy, including an extension of a previously known lower 
bound [22], 

3.1. Collective entropy 

First concepts of information entropies as an uncertainty measure for pointer-based 
simultaneous measurements can be found in Ref. [15], which serves as a foundation for 
the present section. Since information entropies are based on probability distributions, 
we first need to define suitable probability distributions before we can deal with the 
actual entropies. 

The probability of measuring an inferred position X and an inferred momentum V 
is given by the joint probability distribution [13, 46] 


prO* 1 , P; t) = (S(X - X(t))S(V - P(t))) 


( 20 ) 


with the Dirac delta distributions 6(X — X(t)) and S(V — V(t)), which contain the 
inferred observables X(t) and V(t), Eq. (8), respectively. Likewise, the probability of 


measuring either the inferred position X or the inferred momentum V is given by the 
marginal probability distribution of inferred position 


+oo 



(21a) 


— OO 


and inferred momentum 
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+OO 

P <T v {V]t) = JdX pr(X,V]t), (21b) 

— OO 

respectively. 

With the help of the marginal probability distributions, Eq. (21), one can define 
the collective entropy [22] 


S{t) = S x {t) + S v {t), (22) 

which describes the total uncertainty of a simultaneous pointer-based measurement 
process. It consists of the sum of the marginal entropy of inferred position 


+oo 

Sx(t) = - JdX pr x (A; t) In pr^(A; £) 

— OO 

and the marginal entropy of inferred momentum 

+oo 

S P (t) = - J dV pr p(V-,t) lnpr v (V;t). 

— OO 

So far, our considerations have been completely general. 


(23a) 


(23b) 


3.2. Separable initial state 

A more detailed calculation of the collective entropy, Eq. (22), is only possible if we 
choose a more specific initial state £(0) for the measurement configuration. First of all, 
it is a common approach to assume that the bath is initially in thermal equilibrium 
and separable from the system and the pointers |. Furthermore, it seems natural to 
choose initially localized and separable pointer states. These localized pointer states 
are then propagated by the Hamiltonian, Eq. (1), in such a way, Eq. (8), that their 
new location can be read out at end of the interaction process to infer the system 
observables, Eq. (7). A straightforward realization for localized and separable pointer 
states are squeezed vacuum states [47, 3], which have already been used in the classic 
Arthurs and Kelly model [1], Our system to be measured should, on the other hand, 
not be subject to any assumptions, so we describe it by a general density matrix. 

$ For baths which are not initially separable from the system and the pointers, respectively, our 
method of calculating the entropy in Appendix A may fail. However, by choosing an appropriate 
time-dependency of the Hamiltonian, Eq. (1), a smooth switch-on of the bath can be introduced to 
realize a more physically reasonable model; see, e. g., Ref. [10, 23] and references therein. 
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Thus, for our model of open pointer-based simultaneous measurements, we choose 
a separable initial state 


m = £s(0) < 8 > |cri) (ail (8) |er 2 ) (cr 2 | < 8 > £b(0). (24) 

It consists of the initial state £g(0) of the system to be measured, the squeezed vacuum 
states |cri) and |cr 2 ) of the pointers, which can be written as [3] 


(x\o- k ) 



(25) 


in position space with variances of for k e {1, 2}, and the thermal state of the bath 


£b(0) = - exp 




(26) 


with the thermal energy /3 -1 and the normalizing partition function Z. 

The choice of a specific initial state of the bath allows us to determine the statistical 
properties of the stochastic force £(t), Eq. (11). In particular, one has 


Cm = o 


(27) 


and 


\ (£(*i)£ T {h) + £(* 2 )£ T (ti)) = g{ti)g{t 2 )u(ti - t 2 ) (28) 

with the noise kernel [11] 


is{t) 


OO 

J dcucoth 

o 



cos(o;t)I(a;), 


which contains the spectral density [30] 


(29) 


I(cj) = g T m 2 <n 1 5(o;l — o;)m ^g. (30) 

Here we have made use of the Dirac delta distribution <5(a;l —cj) with the identity matrix 
1 and the bath frequency matrix Eq. (12). In brief, the noise kernel describes the 
noisy influence of the bath on the measurement process. Note that we do not further 
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specify the structure of the spectral density § and keep it as a general phenomenological 
expression. Nevertheless, we assume a continuous bath with N —» oo for which the 
definition of a spectral density makes sense in the first place. 

We furthermore remark that for our chosen initial state p(0), Eq. (24), the 
expectation value shift s (t), Eq. (15), vanishes for all times, i. e., 


s (t) = 0. (31) 

This means that the first moments of the inferred observables after a certain interaction 
time t directly correspond to the first moments of the initial system observables, Eq. (14). 
As we have already stated above, this can be considered as a typical behavior. 

3.3. Marginal probability distributions 

The choice of the initial state p(0), Eq. (24), allows us to explicitly calculate the the 
marginal probability distributions, Eq. (21). It seems natural to assume that these 
marginal probability distributions do not directly correspond to the initial position 
distribution (x|ps(0)|ic) or the initial momentum distribution (p|ps(0)|p) of the system, 
but should in addition also incorporate the disturbance from the indirect measurement 
via the pointers as well as the noisy effects of the bath. Since pointers and bath 
are initially in a Gaussian state, Eqs. (25) and (26), their influence on the marginal 
probability distributions can also be expected to have a Gaussian shape. Indeed, the 
calculations shown in Appendix A lead us to the broadened distributions 


pr X ( x >t) 


+oo 



da; (x|^s(0)l a; ) ex P 


(A-x) 2 ' 


(32a) 


and 


pr 


+oo 



dp (p|ds(0) \p) exp 


( V-Pf 

2 5l{t) _ 


(32b) 


respectively with the noise terms S\,(t) and 5j,(t), which represent the diagonal elements 
of the covariance matrix of the noise operators, Eq. (18). As also shown in Appendix 
A, this covariance matrix can be expressed as 

§ A common choice for Eq. (30) is an Ohmic spectral density with I(w) ~ w and a cut-off term for 
high frequencies. See, e.g., Ref. [30] for a more detailed discussion. 
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t t 

+ Jdti j dt 2 g{t 1 )g(t 2 )A{t,t 1 )u(t 1 -t 2 )A T (t,t 2 ). (33) 

o o 

Here we have recalled the time-dependent coupling strength of the bath g(t), Eq. (5), 
the coefficient matrices B (t) and A (t,s) from the dynamics of the inferred observables, 
Eq. (7), the initial value vector J, Eq. (10), and the noise kernel is(t), Eq. (29). Note that 
the non-diagonal terms in Eq. (33) correspond to the correlation term 8xv(t), Eq. (18), 
which is of no further importance in the following. 

Briefly summarized, the marginal probability distributions, Eq. (32), both consist 
of a convolution of the system’s initial probability distributions with a Gaussian noise 
function. This noise function itself is determined by a convolution of a Gaussian pointer 
noise function, Eq. (A. 13), with a Gaussian bath noise function, Eq. (A.22). In other 
words, both the pointers and the bath act as a Gaussian filter [15] through which we 
are forced to look during the measurement process, and which leave us with a distorted 
image of the initial probability distributions of the system. The noise terms 8\{t) and 
8f,{t) give a description for this combined disturbance of the measurement results from 
the pointers and the bath. Since we use Gaussian-shaped initial states for the pointers 
and the bath, Eqs. (25) and (26), which are fully characterized by their second moments, 
the noise terms also contain only second moments: The first term on the right-hand side 
of Eq. (33) describes the pointer-based variance contributions to the noise terms, whereas 
the second term describes the bath-based variance contributions. Both increased 
pointer-based variance contributions and increased bath-based variance contributions 
increase the effective disturbance of the marginal probability distributions, Eq. (32). 
In this sense, the noise terms can also be considered as “extrinsic” or “measurement 
uncertainties” in comparison with the “intrinsic” or “preparation uncertainties” of the 
system given by (x|ps(0)|:r) and (p|(?s(0)|p), respectively. This classification of the 
disturbance is a concept similarly used for variances as uncertainties of pointer-based 
simultaneous measurements [16, 4], 

An optimal measurement configuration for a pointer-based measurement is 
consequently defined by a minimal noise term product. In Appendix B we show that 
the product of the noise terms is bound from below by || 

8 x (t)8 v (t) > -, (34) 

| We remark that we have assumed in Ref. [4] (in Eq. (55b)) that the pointer-based variance product, 
i. e., the product of the diagonal elements of the first term in Eq. (33), is bound from below by 1/4. 
While this assumption always holds true for a closed pointer-based measurement due to Ineq. (34), it 
is in fact not generally valid for an open pointer-based measurement. 
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which defines the best possible accuracy of any measurement apparatus. This statement 
is closely related to the fact that the variance-based uncertainty of a pointer-based 
measurement is bound from below by one, see, e.g., Ref. [1, 31, 20, 4], 

3.4. Lower bound of the collective entropy 

In Ref. [22] we have discussed the collective entropy for closed pointer-based 
simultaneous measurements with pure initial system states and, with the help of 
Ref. [48], have established a lower bound of this collective entropy. Interestingly, 
although we have used an open-pointer based measurement and possibly mixed initial 
system states in the present manuscript, the structure of the marginal probability 
distributions, Eq. (32), which determine the collective entropy, Eq. (22), is similar to the 
structure of the marginal probability distribution for closed pointer-based simultaneous 
measurements. Therefore, with only slight changes, we can adapt the derivation of a 
lower bound of the collective entropy of closed pointer-based simultaneous measurements 
to the collective entropy of open pointer-based simultaneous measurements. In the 
following, we will briefly recapitulate this derivation. 

First of all, we recall a theorem from Ref. [48], which states that the information 
entropy 


of a Fourier convolution 


+OO 


S[f] = - I dx /(xi In fix) 


(35) 


+ OO 

(f* 9 )(x)= J dyf(x-y)g(y) (36) 

— OO 

of two probability distributions f(x) and g(x) is bound from below by 

Slf* 3} > AS[/) + (1 - A)%] - A1 ‘ lA + (1 ~ A) ‘ n(1 - A) (37) 

with an arbitrary weighting parameter A G [0,1]. Equality in Ineq. (37) holds true if 
and only if both f{x) and g(x) are Gaussian distributions with variances a'j and cr|, 
respectively, and the weighting parameter reads 


A 


01 


a s + a l 


(38) 


It is clear that the marginal entropies, Eq. (23), which contain the marginal probability 
distributions, Eq. (32), also represent information entropies of convolutions as given by 




Entropic uncertainty bound for open pointer-based measurements 


15 


the left-hand side of fneq. (37). Therefore, we can apply the lower bound of fneq. (37) 
to each of the marginal entropies in the collective entropy S(t), Eq. (22). In particular, 
we insert the initial position and momentum distributions of the system in place of f(x) 
and their associated Gaussian filter functions in place of g(x). For reasons of simplicity, 
we use the same weighting parameter A for both of these lower bounds. In a next step, 
we can eliminate the dependency on the system to be measured by making use of the 
entropic uncertainty relation [49, 41, 43, 50] 


S[(z|£? s (0)|z)] + S[(p|£s(0)|p>] > 1 + hi7r. (39) 

It is based on the Babenko-Beckner inequality [51, 42], Saturation in Ineq. (39) occurs 
only for pure Gaussian states [52, 53]. For a detailed discussion of this and similar 
entropic uncertainty relations, also see Ref. [45]. The usage of Ineq. (39), which contains 
the density matrix £g(0) of the system to be measured, is the main difference to the 
derivation from Ref. [22], where the system to be measured was limited to pure states 
and a simplified form of Ineq. (39) with pure states had been used. 

Accordingly, we arrive at the lower bound 

S(t) > 1 — A In ^ + (1 — A) In ^ 27r M*)M*) ^ ( 40 ) 

of the collective entropy S(t), Eq. (22). In particular, this lower bound depends on the 
noise terms Sx(t) and S-p(t), Eq. (33). A maximization of the right-hand side of Ineq. (40) 
with respect to A reveals the optimal weighting parameter A = 1/(1 + 2 5x{t)5p(t)) and 
finally leads us to the result 


S(t) > 1 + In 


2vr (5x{t)Sv{t) + 


(41) 


This entropic uncertainty bound for open pointer-based simultaneous measurements is of 
the same form as the entropic uncertainty bound for closed pointer-based simultaneous 
measurements from Ref. [22], Moreover, it is an extension of the more well-known 
entropic uncertainty bound S(t) > 1 + ln(27r) from Ref. [15], to which it can be reduced 
in case of minimal noise terms, Ineq. (34). 

Due to the equality conditions of Ineqs. (37) and (39), equality in Ineq. (41) occurs 
only for initial system states £?s(0), which are minimal uncertainty states in the sense of 
Heisenberg’s uncertainty relation [54], i.e., pure Gaussian states for which the position 
variance a 2 and the momentum variance a 2 obey o' 2 a 2 = 1/4; and which additionally 
fulfill cr 2 = 5x(t)/(25-p(t)) for a given set of noise terms [21]. Such initial system states 
can be understood as “minimal entropy states” [22]. Note that mixed system states 
with Gaussian density matrices [55] are not sufficient for equality. 

Summarized, the collective entropy of an open pointer-based measurement, 
Eq. (22), behaves exactly like the collective entropy of a closed pointer-based 
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measurement. Particularly, the collective entropy is bound from below by the sharp 
entropic uncertainty bound given by Ineq. (41). The only effect of the environmental 
heat bath is to modify the collective entropy by modifying the noise terms, Eq. (33). An 
optimal measurement accuracy can be achieved for minimal noise terms, Ineq. (34), and 
requires a pure Gaussian system state with a specific variance which leads to equality 
in Ineq. (41). 

4. Conclusion 

We have shown that it is possible to determine a formal expression for the collective 
entropy of a general open pointer-based simultaneous measurement. This collective 
entropy has a very intuitive structure from which the noisy influence of the pointers 
and the bath on the measurement result can be understood. Moreover, this structure 
allows us to show that the collective entropy has a sharp lower bound, which can only 
be reached for specific pure Gaussian system states. In particular, our results are valid 
for any bilinear interaction between the system and the pointers and any initially mixed 
system state and thus extend various previous results on this topic. 

Several simplifications have been made to perform our calculations. First of all, the 
Hamiltonian only includes bilinear terms. However, we expect that terms of higher order 
would only lead to relatively small correction terms without fundamentally changing 
our results. Second, ideal single variable measurements of the pointer observables have 
to be performed in order to realize the measurement procedure, which is an inherent 
conceptional weakness of pointer-based simultaneous measurements. Yet we hope that 
with the help of the environmental heat bath, it might be possible to replace this rather 
theoretical measurement process by a more natural decoherence process [56]. Finally, we 
have used a separable initial state with squeezed states as the initial pointer states and 
a thermal state as the initial bath state. Although this approach is a clear limitation of 
our results, we think that our choice is reasonable. It would nevertheless be interesting 
to discuss different initial states and specifically examine the influence of entanglement 
and preparation energy on the collective entropy. 
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Appendix A. Marginal probability distributions 

In this appendix section we briefly describe how to calculate the marginal probability 
distributions, Eq. (32), from the definition, Eqs. (20) and (21), with the help of the 
chosen initial state, Eq. (24). Our calculations are mainly based on characteristic 
functions and their Fourier transforms. 
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Inserting the formal dynamics of the inferred observables, Eq. (7), 
probability distribution, Eq. (20), leads to 


into the joint 


pr(A’,'P; t) = tr 


x 8 


X - X s (0) - (1, O)B(t) j - (1, 0)(A x £)(t) 
V - Ps(0) - (0, l)B(t) J - (0,1)(A X i)(t) 


(A.l) 


The Dirac delta distributions in Eq. (A.l) can be expressed as Fourier transforms of 
unity, i. e., 8{x) = T {1} (x). Here we use the notation 


•H/K*) 


i 

27T 


+00 

J da exp[i ax] /(a) 


— OO 


(A.2) 


for the Fourier transform of an arbitrary function f(x), and an analogous notation for 
Fourier transforms of functions of two variables. Thus, the separability of the initial 
state £?(0), Eq. (24), allows to rewrite Eq. (A.l) as 


pr(A, V, t) = F {F S F P F B } (A, V) = [T {F s } * F {F P } * 7 {F B }) (A, V) (A.3) 
with the characteristic function of the system 

F s (a 1 ,a 2 ) = tr |^ s (0)exp[-i{aiA s (0)+a 2 P s (0)}]| , (A.4) 

the characteristic function of the pointers 

F P (ai, a 2 ) = (cti| (g) (<j 2 \ exp[— i{(ai, a 2 )B(t)J}] |ai) ® |cr 2 ), (A.5) 

and the characteristic function of the bath 


F B (cti,a 2 ) = tr |^ B (0)exp[-i{(ai,a 2 )(Ax|)(t)}]| . (A.6) 

The symbols x and * are defined in Eqs. (13) and (36), respectively. 

In the following, we calculate the Fourier transforms of the characteristic functions, 
Eqs. (A.4) to (A.6), one after another with the help of Wigner functions [3]. This allows 
us to explicitly perform the convolution in Eq. (A.3). By straightforward integration 
of the resulting Gaussian expressions we finally arrive at the marginal probability 
distributions, Eq. (32). 
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Appendix A. 1. System 

First of all, the Fourier transform of the characteristic function of the system, Eq. (A.4), 
corresponds to the respective Wigner function 


W s (x,p)=F{F s }(x,p) (A.7) 

of the system state. We can also use this connection between characteristic function 
and Wigner function to determine the Fourier transforms of the other two characteristic 
functions, Eqs. (A.5) and (A.6), in the following. 

Appendix A.2. Pointers 

The Wigner function W a of a squeezed vacuum state |cr) with position variance a 2 can 
be written as 


W a (x,p) = G ^-^,-2u 2 ,0, , (A.8) 

where we have introduced the general Gaussian function 

G (a, b, c, n; x,p) = n exp [ax 2 + bp 2 + exp] . (A.9) 

Equation (A.8) corresponds to the Fourier transform of the generic characteristic 
function 


F^(aq,a 2 ) = (ct| exp[-i(aqx + a 2 p)]|er) (A.10) 

of squeezed vacuum states as we have used them for the pointer states, Eq. (25), where 
x and p stand for the position and momentum observables, respectively, in the Hilbert 
space of | cr). Looking at this relation the other way round yields 

K( X ,P) = f- 1 {W a }(x,p) = G ^-y,-^,0,l;x,p^ , (A.11) 

where J 7-1 {/} ( x ) denotes the inverse Fourier transform of an arbitrary function ffx) 
with F {F^ 1 {/}} (x) = f{x), Eq. (A.2). Thus, the characteristic function of the 
pointers, Eq. (A.5), can be written as 


F P (a 1 ,a 2 ) = F^((ai,a 2 )B(f)(l,0,0,0) T , (aq, a 2 )B(t)(0, 0,1,0) T ) 

x ^((cq, a 2 )B(t)(0,1, 0, 0) T , (cq, a 2 )B(t)(0, 0, 0,1) T ) (A.12) 
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T{F?}(x,p) = G 


ap(t) bp(t ) c P (t) 1 
dp(t) 1 dp{t) ’ dp{t) ’ dp(t) 


x,p 


with the coefficients 


(A.13) 


and 


— 2bp{t) cp(t) \ 

cp(t) — 2ap(t) J 


B{t)\B T (t) 


(A.14a) 


dp(t) = Aap(t)bp(t) — Cp(t), 


where 


V = \ ({JJ T > + (Jj t > t ) 


K 

0 

0 

0 ^ 

0 


0 

0 

0 

0 

1 

4-? 

0 

1° 

0 

0 

4 k) 


(A. 14b) 


(A.15) 


denotes the symmetrized pointer covariance matrix, which directly follows from the 
definition of the initial value vector J, Eq. (10), and the initial pointer states, Eq. (25). 
In conclusion, Eq. (A. 13) can be calculated solely from the coefficient matrix B(t), 
Eq. (7), and the initial pointer position variances o\ and cr|, Eq. (25). 


Appendix A.3. Bath 

The Wigner function of the thermal bath state £>b( 0), Eq. (26), reads 


W B (q,k) = det 

x exp 


tanh ( ) 7r 1 


u>l3 


m 2 tanh - ojrim q 


1 rp _/ LO j3 . _ - _ . , 

k m 2 tanh - oj m 2 k 


(A.16) 


Analogously to the calculation of the characteristic function of the pointers, the generic 
characteristic bath function 
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^b(“i>“ 2 ) = tr |^b( 0) exp[— i(a^q + a^k)]| (A.17) 

of baths in thermal equilibrium, which is connected to the characteristic bath function, 
Eq. (A.6), via 

F b (oi i, a 2 ) = F g ((A * u q ) T (t) (ai, a 2 ) T , (A * u k ) T (t) (c>q, a 2 ) T ) , (A.18) 

can be written as 


F'(q,k)=^- 1 {fE B }(q,k) 
= exp 


1 rp 1 f CO f3 \ -l _ 1^ 

-q m 2 coth - co m 2 q 

2 1 H 


1, T I , (U/3 \ 1 

-k m 2 coth - wrrrk 

4 2 


Here we have made use of the abbreviations 


u q (t) = —g T (f)m 2 cos(u>f)m 


and 


(A.19) 


(A.20) 


u k (i) = — g T (f)m 2 sin(u;t)u; 1 m 2 . 
Explicitly performing the Fourier transform of Eq. (A.18) leads to 


(A.21) 


E{F B }(x,p)=g 


Qb(£) b B (t) c B (t ) 1 

d B (t) 1 d B (t) ’ d B (t) 1 2ir^/d B (t) 


\x,p 


with the coefficients 


(A.22) 


-2 b B (t) c B (t ) \ , w,w A , Tu \ 
c B (t) -2« B ( t )h (A * v)W(A * v) (t) 


(A.23a) 


and 


d B (t) = 4a B (t)b B (t) - c B (t), 


(A.23b) 
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v(ti)v r (t 2 ) = g{ti)g(t 2 )is(ti - t 2 ) (A.24) 

with the noise kernel i/(i), Eq. (29). 

As a result, Eq. (A.22) is determined by the coefficient matrix A(t,s), Eq. (7), 
and the noise kernel i/(£), Eq. (29). Note that for a turned off bath (i. e., g (t) = 0 
for all t > 0), one has d B (f) = 0 and thus Eq. (A.22) is not well-defined. However, 
Eq. (A.22) can in this case be understood as a Dirac delta distribution so that the 
following considerations are still applicable with a p (t) = 0, b p (t) = 0 and c B (£) = 0. 

Appendix A.f- Coefficients 

In a next step, we show how the pointer coefficients a P (t), b P (t ) and cp(f), Eq. (A. 14a), 
and the bath coefficients a B (t), & B (f), and c B (f), Eq. (A.23a), are related to the noise 
terms 8^{t) and 8f>{t) and the correlation term 5xp(t), Eq. (18). For this purpose, we 
recall the noise operators Nx{t ) and N-p(t), Eq. (16). According to Eqs. (17) and (31), 
we can simplify their covariance matrix, Eq. (18), to 


( %(t) SMt)\ = /( N%(t) 

V-ww sm) \\NMt) Np t ))/■ 1 - ; 

Here we have also recalled the symmetrized noise operator Nxvif)-, Ed- (19). A 
straightforward calculation using the definition of the noise operators, Eq. (16), shows 
that 


(<W(f) ^pt)) = B WVB T («) + (A*v)(«)(A*v) T (i) (A.26) 

with the abbreviations introduced in Eqs. (A. 15) and (A.24). Here we have also made 
use of the symmetry relation [11] 


m yem = (meihf (A.27) 

for the stochastic force £(t), Eq. (11). A comparison of Eq. (A.26) with Eqs. (A.14a) 
and (A.23a) immediately reveals 

( &x(t) Sxr{t)\ = f-2(bp{t) + b B {t)) c P (t) + c B (t) \ 
y>xr{t) 8p(t) J y cp(t) + c B (t) —2(a P (t) + a B (t)) J 


which relates the noise and correlation terms on the left-hand side with the pointer and 
bath coefficients on the right-hand side. 
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Appendix A.5. Probability distributions 

At this point, we can collect onr previous results. By inserting Eqs. (A.7), (A. 13) 
and (A.22) into Eq. (A.3), we arrive at the final expression for the joint probability 
distribution 


p(A.29) 
with the coefficients 


A 2 (t) ~ d ® 

2 (ap(t) + a B (t)) ’ 

(A. 30a) 

. 9 . . ~d(t) 

p{ ) ~ Hb P (t)+b B (t))' 

(A.30b) 

- cp(t) + c B (t) 

l(t) ~ ax, • 

(A.30c) 


and 


d{t) = 4 (a P (t) + a B (t))(b P (t) + b B (t)) - (c P (t) + c B (t)) 2 . (A.30d) 

Finally, the marginal probability distributions, Eq. (32), follow from Eq. (A.29) by 
straightforward integration as defined in Eq. (21) when we make use of the marginals 


and 


+oo 

J dpW s (x,p) = (x|£ s (0)|x) 

— OO 


(A.31a) 


+oo 

J d xW s (x,p) = {p\ps(0)\p) (A.31b) 

— OO 

of the initial Wigner function of the system Ws(x,p), Eq. (A.7), where (x|£s(0)|a;) 
and (p|^s(0)|p) stand for the initial position distribution and the initial momentum 
distribution of the system state, respectively. In particular, the resulting marginal 
probability distributions, Eq. (32), contain the noise terms d 2 x (t) and Eq. (A.28). 
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Appendix B. Minimal noise term product 

We show in this appendix section that the product of the noise terms 8\{t) and S^(t), 
Eq. (18), obeys Ineq. (34). Our proof is closely related to the considerations in Ref. [31]. 
Due to the definition of the noise terms, Eq. (18), we can make use of Robertson’s 
uncertainty relation [57], to establish the lower bound 


Sx{t)S-p(t) > - 


N x {t),N v (t) 


) 


which contains the commutator 


(B.l) 


N x (t),N v (t) = X{t),V{t) - *(t),P s (0) - A S (0 ),V{t) + A s (0),P s (0) 


(B.2) 


of the noise operators N x (t ) and N-p(t), Eq. (16). The first commutator on the right- 
hand side of Eq. (B.2) vanishes due to Eq. (9). The second and third commutators 
read 


X(t),P s (Q) 


*s(0 ),V(t) 


A's(0).P s (0) , 


(B.3) 


since only the first term on the right-hand side of Eq. (7) acts in the initial system’s 
Hilbert space. As a result, we have 


N x (t),N v (t) 


A s (0),P s (0) 


(B.4) 


Inserting Eq. (B.4) into Ineq. (B.l) finally leads to Ineq. (34). 

We remark that it would also be possible to use Schrodinger’s uncertainty relation 
[58] instead of Robertson’s uncertainty relation, Ineq. (B.l). This approach would 
additionally incorporate the correlation term S x -p(t), Eq. (18), on the right-hand side of 
Ineq. (34), i.e., 


$(*)$(«) > s'Ut) + \- 


(B.5) 


We do not discuss this extension to Ineq. (34) in more detail. It could, however, serve 
as an interesting point of origin for further considerations. 
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